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We use supercomputer simulations to show that inter-atomic interactions can strongly affect the 
phase evolution of Bose-Einstein condensates that are diffracted from atom chips, thereby explaining 
recent experiments. Interactions broaden and depopulate non-central diffraction orders and so, 
counter-intuitively, they actually reduce the spatial width of the diffracted cloud. This means, 
more generally, that the experiments cannot be fully described by non-interacting theories and, 
consequently, interferometers that use many interacting atoms to enhance sensitivity require complex 
calibration. We identify device geometries required to approximate non-interacting behavior. 
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I. INTRODUCTION 

The ability to trap Bose-Einstein condensates (BECs) 
close to a surface using microchip technology provides a 
highly versatile way to manipulate and study the quan- 
tum wave nature of cold atoms [H, 0, i, i, S, H, B i, 
[Tol . [III . d, [3, ■ Consequently, atom chips form a 
clean and flexible environment for matter wave interfer- 
ometry [l^, [l^ [13, [H, [l^ . Recent experiments P, [1| cre- 
ated an atom interferometer by reflecting a BEC from a 
magnetic grating fabricated on such a chip. Although ev- 
idence of successful diffraction was obtained, sharp, well- 
separated, diffraction orders were not observed. Further- 
more, to extract phase information from the BEC, the 
experiments relied on two assumptions. Firstly, inter- 
atomic interactions do not change the BECs phase in 
any way. Secondly, the reflection process is short enough 
to modulate only the phase of the BEC, but not its den- 
sity. It is unclear to what extent these assumptions are 
valid and how they affect the accuracy of phase data in- 
ferred from the experiments. 

In this paper, we show that inter-atomic interac- 
tions strongly affect the momentum distribution of the 
diffracted atom cloud, and hence its phase evolution, by 
broadening the diffraction orders and depopulating the 
higher ones. In addition, the reflection process is long 
enough to modulate strongly both the phase and den- 
sity of the condensate. When these effects are included, 
our calculations give good quantitative agreement with 
the diffraction patterns observed in the experiments. In- 
terferometers based on interacting par ticles could make 
sensitive measurement devices [idl. llSl. [20| . but our re- 
sults show that they respond in a highly complex way to 
spatially-varying phase patterns. To recover the essen- 
tially non-interacting behavior required for simple ana- 
lytical models to make accurate a priori predictions of an 
interferometer's characteristics, we find that the kinetic 
energy of the BECs separating diffraction orders must 
exceed the peak inter-atomic interaction energy. 

The diffraction process has two distinct stages: the im- 
printing stage in which the BEC reflects from the grat- 




FIG. 1: (Color) (a) Schematic cross-section (black shape) in 
the a; = plane (axes inset) through the microchip and wires, 
which are separated by d/2. Blue shape shows the BEC in 
its initial position centered at y = 0. Left (right) dotted 
parabola shows the harmonic trap potential before (after) 
displacement Ay. (b) Curved surface: real part of the po- 
tential energy VT{y,z) (defined in the text, axes inset) after 
trap displacement. The displaced trap accelerates the BEC 
(blue) towards the microchip and wire grating on the right 
(axes inset). Vertical arrows on the wires indicate direction 
of current. Curved black arrow indicates the BECs motion 
path. The offset magnetic field is not shown for clarity. 

ing, and the expansion stage in which the atom cloud is 
released from its trap [l|, |3| . The imprinting process en- 
codes a diffraction pattern in momentum space, which, 
after time-of-flight expansion, manifests itself in a real- 
space diffraction pattern. Our goal is to understand the 
role of interactions and excitations at each stage. We 
therefore perform separate simulations for the reflection- 
imprinting process and the subsequent expansion. 



II. SYSTEM AND METHODOLOGY 

Figure [T] shows a schematic diagram of the system used 
to imprint the BEC. In the experiments and our cal- 
culations, a BEC of 1.2x10^ ®^Rb atoms is prepared in a 
harmonic trap with frequencies uJz = 27rxl6 rad and 
LOr = 27rx76 rad s^^ in the longitudinal (z) and radial 
{x-y plane) directions respectively. The peak atom den- 



sity (no « 10^" m"'^) is initially at x = y = z = 0. An 
array of current carrying wires is microfabricated on a Si 
surface that lies in the y — 30 /ini plane [Fig. [Ha)]. Ad- 
jacent wires are separated by 2 /im [Fig. [Ija)] and carry 
a current of 0.2 mA in alternating directions [Fig. [Hb)], 
giving a grating period d = 4 /im. An offset field of 2 G in 
the z-direction combines with the field from the grating 
wires to create a potential profile that oscillates strongly 
between peaks and troughs, thus forming a diffraction 
grating [Fig.[T]^b)] within ~5 fim of the surface (2ll |. 

At time t = 0, we displace the harmonic trap through 
a distance Ay = 11-15 /im in the y-direction and hence 
accelerate the BEG towards the grating along its short 
radial axis [Fig.[lja)]. Since Ay < 30 /im, the BEG must 
rise up the other side of the trap and decelerate before 
making contact with the grating potential. Impact ve- 
locities, Vy (determined by Ay), are in the range - 
4 mm s~^. The imprinting stage lasts 12 ms, by which 
time the condensate has approximately returned to its 
original position. 

We determine the dynamics of the BEG in 3D by solv- 
ing the Gross-Pitaevskii equation (GPE) 

- + Vt-9 + l^-l-^-ift— , (1) 

Zm m ot 

where m is the mass of a single ®^Rb atom, a = 5.4 nm 
is the s-wave scattering length, is the Laplacian, 
Vriy, z) is the potential energy produced by the magnetic 
trap, grating, and chip surface [Fig. [TJb)], and 4'(a;,y, z) 
is the wave function at time t, normalized so that l^"!^ 
is the number of atoms per unit volume. We solve the 
equation using the Grank-Nicolson ^2] method, and also 
use the Truncated Wigner 23] approach to investigate 
the effects of quantum fluctuations. Since the potential 
landscape contains very sharp features [Fig. [TJb)], solv- 
ing Eq. [T] requires parallel programming implemented on 
a supercomputer to achieve the required resolution [2^ . 

III. RESULTS AND DISCUSSION 

A. The imprinting process 

Figure [5] shows density slices of the BEG in the x — 
plane as it reflects from the wire array following a trap 
displacement Ay = 13 fun. Figure [2I^a) shows the initial 
state at i = ms. Figure [2jb) shows the BEG just after 
contact with the lattice potential aX t = 5.6 ms and Fig. 
mjc) shows the BEG as it begins to move away from the 
lattice at t = 7.6 ms. The atoms interact strongly with 
the lattice for a time interval Af ~ 1 - 2 ms, which de- 
pends on Vy. By the end of this time interval, a density 
imprint and a clear phase imprint [2^ have formed [Fig. 
^c) plus inset]. At t = 12.0 ms, just before the trap 
is switched off, the density imprint is still evident [Fig. 
[DJd)]. For all trap displacements, our calculations reveal 
phase and density imprints that vary approximately si- 
nusoidally across the reflected BEG. 
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FIG. 2: (Color online) Density slices (dark high) in the a; = 
plane [axes in (b)] of the diffracting BEC at t = (a) 0.0 ms, 
(b) 5.6 ms, (c) 7.6 ms and (d) 12.0 ms. Bar in (a) shows scale. 
Black shapes (not to scale) indicate the position of the atom 
chip and grating. Inset in (c) shows phase imprint along z 
(vertical) at t = 7.6 ms. 

Previous work predicts that BEG reflection processes 
can produce vortices in certain parameter regimes (2^ . 
[26l.[27||. Here, however, no vortices form for any trap dis- 
placements because At is always shorter than the time 
required for sound waves to propagate to the edge of the 
condensate and thus trigger the formation of excitations 
[2^ . In addition, quantum fluctuations have been shown 
to produce decoherence and scattering halos in reflected 
BEGs. For the parameters considered here, though, they 
make no significant difference to the BEG's dynamics be- 
cause the velocities involved are too low to generate sig- 
nificant incoherent scattering 28]. The bare GPE there- 
fore captures the key physics and so we use it to produce 
all the results presented here. 

B. Broadening and depopulation of non-central 
diffraction orders 

We now perform a detailed analysis of the role of in- 
teractions in the diffraction process. To isolate their ef- 
fects, we consider a simple model, which does not de- 
pend on the details of the imprinting process, only on the 
BEG's mean-field response to a sinusoidal phase imprint 
4> = S s'm{27rz/d), like that produced by the grating. We 
use S — lA for all our expansion simulations as this is a 
typical value deduced from experiment P, 0] ■ 

Figure [3] shows how the density profile along z evolves 
after phase imprinting a non-interacting atom cloud (left- 
hand column) and an interacting BEG (central column). 
The initial density profiles without and with interactions 
[Fig.[2Ia) and[31[b) respectively] are taken to be identical 
to simplify the comparison of their subsequent evolution. 
The third column in Fig. [3] shows fc-space momentum 
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FIG. 3: Left (center) column: density profiles along the z-axis 
for a non-interacting atom cloud (interacting BEG) after a si- 
nusoidal phase imprint. Right column: Fourier power spectra 
calculated for the interacting BEG. The first, second, third 
and fourth rows down show, respectively, profiles at ms, 
0.2 ms, 1 ms, and 8 ms after the imprint. Dashed curve in 
(1) shows the momentum distribution for diffraction from a 
grating with half the original period (see text). 



distributions for the interacting BEC, obtained by inte- 
grating the wavefunction along the x- and y-directions 
and calculating the one-dimensional Fourier power spec- 
trum, F, along the z-direction. The initial momentum 
distribution has sharp peaks [Fig. [3{c), arrowed] at the 
0th, 1st and 2nd diffraction orders. 

The second row of Fig. [3] shows the density profiles 
and momentum distribution 0.2 ms after the phase im- 
printing. The density profiles for the non-interacting 
[Fig. EJd)] and interacting case [Fig. [3{e)] now contain 
sharp modulations, which are interference fringes be- 
tween emerging diffraction orders. Both density profiles 
look very similar and the momentum distribution for the 
interacting case [Fig. [Hf)] has barely changed from its 
initial form [Fig. EIc)], indicating that the interactions 
have a negligible effect on this timescale. This is because 
0.2 ms is comparable with the minimum correlation time, 
tc = 0.14 ms [29|, which means that collective interaction 
effects have had insufficient time to develop. 

Figure[3]reveals that we cannot, as in previous analyses 
P, , assume that there is no change in the atom density 
profile during the imprinting stage. Since Ai ~ 1 - 2 ms. 



it is not meaningful to speak of a pure phase imprint in 
the experiments [ll because a combination of phase and 
density modulations inevitably occurs on this timescale. 
It may, however, be possible to quantify a heuristic "over- 
all imprint" strength. To obtain a pure phase imprint 
requires At < 0.1 ms. This could be achieved by flashing 
the grating on briefly, by turning the current on and off 
rapidly, as opposed to having it on permanently. 

The third row of Fig. [3] shows the density and momen- 
tum profiles 1 ms after imprinting. Differences between 
the non-interacting and interacting cases now become ap- 
parent because the system has evolved for several corre- 
lation times. In the non-interacting profile [Fig. [3Ug)], 
the density minima extend all the way to zero whereas 
they generally do not in the interacting case [Fig. ^h)] 
[30t . The amplitude of the density modulation is also 
much lower in the interacting case [note the change of 
vertical scale between Figs. |3Ig) and[3Ih)]. The momen- 
tum distribution of the BEC [Fig. [3)^1)] reveals dramatic 
depopulation of the 1st and 2nd order diffraction peaks 
[arrowed]. Counter-intuitively, the repulsive interactions 
therefore reduce the width of the BEC in momentum 
space. This demonstrates that the experimental data 
cannot be fully analysed using a non-interacting theory. 
Such theories assume that the order populations do not 
change during diffraction; this is clearly violated in the 
interacting case. 

Depopulation occurs in the interacting case because 
mean-field repulsion pushes atoms away from high den- 
sity regions into the gaps between the peaks. This pre- 
vents the density minima from falling close to zero [Fig. 
|3l[h)] and broadens the modulations, thereby reducing 
the number of atoms in the higher momentum diffraction 
orders [Fig. [2{i)]. In the non-interacting case, there is no 
mean-field resistance to sharp modulation of the density 
profile and the minima fall close to zero [Fig. [SJg)]. In- 
terestingly, interactions do not significantly change the 
positions of the diffraction orders along the kz axis [com- 
pare Figs, ^c) and[3^i)], and hence the speeds at which 
distinct diffraction orders separate are largely unaffected. 

The fourth row of Fig.[3]shows the density and momen- 
tum profiles at 8 ms after imprinting [solid curves] . The 
non-interacting density profile [Fig.[3jj)] has retained its 
sharp modulations, which extend down to zero. In the 
interacting case [Fig. [3fk)], however, the density mod- 
ulations have become even broader and lower in ampli- 
tude. At this point in the expansion process, the den- 
sity of the BEC has become so low that interactions no 
longer play a significant role, meaning that the diffraction 
peaks subsequently separate at a speed determined sim- 
ply by their positions in momentum space. Consequently, 
the corresponding momentum distribution [solid curve in 
Fig. EI^l)] has the same shape as the far-field spatial den- 
sity profile. It reveals further depopulation of the higher 
diffraction orders and broadening of all the peaks. This 
shows that the experiments could never observe spatially 
well-separated diffraction peaks, even for very long ex- 
pansion times (> 100 ms). Quantum fluctuations were 
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FIG. 4: (Color) (a)-(c) Density profiles along the z-axis after 
a sinusoidal phase imprint and 26 ms expansion time for the 
non-interacting (a) and interacting (b) clouds, (c) Curve ob- 
tained after passing the trace in (b) through a low resolution 
filter, (d) Density profile of the interacting cloud in the y-z 
plane (axes inset), obtained by integrating |^'(a;, z)^ along 
X after a sinusoidal imprint, 26 ms expansion and application 
of a low pass filter to simulate experimental image resolution. 
Horizontal bar shows scale. Density profile in (c) is taken 
along dotted line in (d). 



Illjb) through a low-pass filter with the same resolution 
as the camera used to image the BECs in experiment 
(pixel size ~ 5 /im) . Figure [3Id) shows the correspond- 
ing (filtered) 2D density profile in the y-z plane after 
integrating |^'(a;, y, z)p along x. We now obtain good 
quantitative agreement with the experimental data [com- 
pare Fig. m^d) with Fig. 2(a), right-hand panel, of Ref. 
[3]]. Our simulated density profile has approximately the 
same shape and width as the experimental profile, and 
both contain three small peaks [dotted arrows in Fig. 
Ufc)], corresponding to aliased interference fringes be- 
tween separating diffraction orders. As in Fig. [31 these 
peaks are dramatically weakened compared with the non- 
interacting case where nodes are observed [Fig.[3]^g)]. We 
also reproduce the "shoulders" [labelled "1" in Fig. SJc)] 
observed in experiment, which are remnants of the 1st 
diffraction orders. Since our zero-temperature mean-field 
theory agrees well with experiment, we infer that possi- 
ble deleterious effects due to heating of the atom cloud 
during diffraction are not significant. 



D. Achieving a non-interacting regime 



again found to make no significant difference to these 
simulations, with coherence losses always '--^ 1%. 



C. Comparison with experiment 

To compare our theory with the diffracted BEG den- 
sity profiles seen in the experiments , we now perform 
longer expansion simulations, which approximate the ex- 
perimental procedure. After imprinting both phase and 
density profiles similar to those shown in Fig. [H^c), the 
BEG is held in the harmonic trap for 6 ms and then the 
trap is turned off, allowing the BEG to expand freely for 
20 ms. Figures m^a) and[3]Jb) show density profiles along 
the z-axis at the end of this process for a non-interacting 
matter wave (a) and a BEG (b) . For the non- interacting 
atom cloud [Fig.[4l^a)], the diffraction peaks separate and 
a typical far-field diffraction pattern appears. The 0th, 
1st and 2nd diffraction orders (arrowed) are all signifi- 
cantly populated. For the interacting BEG [Fig.lH^b)], 20 
ms expansion is insufficient to reach the far-field regime. 
Hence the density profile has rapid modulations along the 
z-direction, which are interference fringes between differ- 
ent diffraction orders. As explained in our discussion of 
Fig. [3l^h) , the amplitude of the fringes has been reduced 
by the inter-atomic interactions. Figure ID^b) also shows 
that the spread of the cloud is significantly less than 
the non- interacting BEG. This emphasises the counter- 
intuitive effect of repulsive interactions in a diffraction 
system and suggests that attempts to infer phase data 
from atom interferometers based on BEGs may give rise 
to major systematic errors if interactions are not properly 
accounted for. 

Figuredl^c) shows the effect of passing the curve in Fig. 



We now consider the parameter regimes required to 
approximate non-interacting behavior. We find that we 
approach this regime gradually as the kinetic energy asso- 
ciated with the first diffraction order, Ed = iir^h^ /md? 
[3l| . becomes comparable to, or greater than, the BEGs 
peak interaction energy Ei = 47rh?ano/m. This condi- 
tion can be satisfied by halving the original grating pe- 
riod to d = 2 /im. In this case, after 8 ms expansion, the 
diffraction peaks [dashed curve in Fig. EI^l)] are very sim- 
ilar to the non-interacting case [Fig. [Sljc)] and reveal far 
less depopulation and broadening than when d = A /im 
[solid curve in Fig. ^\)]- We obtain similar results by 
keeping d — 4 /im but reducing the BEG's density so 
that the condition Ejy > Ej is satisfied. An alterna- 
tive approach is to increase radial confinement so that 
the BEG expansion occurs on a faster timescale than the 
depopulation of the diffraction orders. 



IV. CONCLUSION 

In summary, we have studied how inter-atomic interac- 
tions affect the dynamics of BEGs diffracting from mag- 
netic gratings on atom chips. Analysis of our numerical 
calculations reveals that if Ej > Ed, as in recent ex- 
periments P, inter-atomic interactions significantly 
broaden the diffraction orders and depopulate the higher 
ones. This means that the BEG's phase evolution cannot 
be described accurately by simple non-interacting theo- 
ries, which omit key physical processes and therefore lead 
to qualitatively incorrect predictions. Atom chips like the 
grating considered here, which imprint well-defined pat- 
terns in the phase and density profiles of a BEG, are com- 
monly used as atom interferometers. Our results show 
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